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The geometric Levy model (GLM) is a natural generalisation of the geometric 
Brownian motion model (GBM) used in the derivation of the Black-Scholes formula. 
The theory of such models simplifies considerably if one takes a pricing kernel ap- 
proach. In one dimension, once the underlying Levy process has been specified, the 
GLM has four parameters: the initial price, the interest rate, the volatility, and the 
risk aversion. The pricing kernel is the product of a discount factor and a risk aver- 
sion martingale. For GBM, the risk aversion parameter is the market price of risk. 
For a GLM, this interpretation is not valid: the excess rate of return is a nonlinear 
function of the volatility and the risk aversion. It is shown that for positive volatility 
and risk aversion the excess rate of return above the interest rate is positive, and is 
increasing with respect to these variables. In the case of foreign exchange, Siegel's 
paradox implies that one can construct foreign exchange models for which the excess 
rate of return is positive both for the exchange rate and the inverse exchange rate. 
This condition is shown to hold for any geometric Levy model for foreign exchange 
in which volatility exceeds risk aversion. 
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I. INVESTMENT-GRADE ASSETS AND EXCESS RATE OF RETURN 

The goal of this paper is to understand better the nature of the "risk premium" associated 
with jumps in asset prices. The idea is to work in a rather general setting, without being 
tied too much to any particular model. For mathematical simplicity, we shall assume that 
the dynamics of asset prices are driven by Levy processes. This already encompasses a large 
class of models — it includes, above all, all the Brownian motion based models — so we need 
not worry that we are being too restrictive. We aim to clarify the relation between risk, 
risk aversion, and the excess rate of return (above the interest rate) offered by risky assets 
in such a context. With this end in mind, let us recall the setup in the geometric Brownian 
motion (GBM) model. The GBM model is very simple, but it captures a number of the 
main features of the relation between risk, risk aversion, and the excess rate of return. We 
shall adopt a pricing kernel approach, which turns out to be particularly useful in the Levy 
case since it allows one to distinguish clearly between pricing issues and hedging issues. See, 
for example, Cochrane (2005) for an overview of the application of pricing kernel methods. 

In the one-dimensional case, we have a Brownian motion {Wt}t>o 011 a probability space 
(Q, T, P), and the associated augmented filtration {Ft}t>a- Here P represents the physical 
measure, and {Ft} is taken to be the market filtration. The model consists of: (a) a pricing 
kernel; and (b) a collection of one or more "investment-grade" assets. For simplicity, we 
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assume for the time being that the assets under consideration pay no dividends over the 
time horizon considered. We relax this assumption in Section VI. 

The idea of an investment-grade asset is that it should offer a positive excess rate of return 
above the interest rate. There are respectable assets that do not have this property — such 
assets are typically held alongside investment-grade assets for hedging. One can check, for 
example, in the context of the Black-Scholes model, that the price process of a put option 
has a drift that is less than the interest rate. The pricing kernel in the GBM model is 

n t = e -rt e -AWi-^» tj (1) 

where r > is the interest rate, and A > is the risk aversion factor, both assumed to be 
constant. For a typical investment-grade asset we then have 

S t = S e {r+Xa)t e aWt -^ 2 \ (2) 

where a > is the volatility. The term Xa is called the "risk premium" or "excess rate of 
return" , and is clearly positive under the assumptions made. 

We observe that the risk premium is increasing with respect to both the volatility and 
the risk aversion. Since Xa is linear in each factor, we call A the "market price of risk" in 
the GBM model. It should be evident, however, that there is no a priori reason why the 
excess rate of return should be bilinear. Indeed, we shall demonstrate that in a general Levy 
model the excess rate of return is a nonlinear function of A and a. The reason that the 
pricing kernel is such a useful concept in finance is that market equilibrium and the absence 
of arbitrage are both built into the idea that the product of the pricing kernel with the price 
of any asset paying no dividend is a martingale. In the GBM case, for example, we have 

n t S t = S e( CT - A)m ^ (CT - A)2 '. (3) 

We shall use this property of the pricing kernel to establish the general form of an arbitrage- 
free Levy-driven asset-pricing model. In Section [XT] we look at one-dimensional geometric 
Levy models with constant volatility and risk aversion, and in Section II III we establish 
the positivity and monotonicity of the risk premium for such models. In Section IIVI we 
consider models for foreign exchange, and establish conditions sufficient to ensure that both 
the exchange rate and the inverse exchange rate have a positive risk premium. Examples of 
geometric Levy models are studied in detail in Section |Vj where we note the fact that, unlike 
the GBM case, option prices in general depend on the risk aversion level; then in Section IVT1 
we extend the model to include dividends. In Section IVHI we consider models in which the 
market is driven by a vector of Levy processes, and in which the volatility and risk aversion 
are predictable processes. In both situations we establish conditions sufficient to ensure the 
positivity of the excess rate of return. In conclusion we comment on the advantage of the 
use of pricing kernel methods in the theory of Levy models, and how this allows us to unify 
earlier work on the subject, leading to a coherent framework for asset pricing. 



II. PRICING KERNEL APPROACH TO LEVY MODELS FOR ASSET PRICES 

Let us construct a family of Levy models in the spirit of the GBM model. We shall call 
these geometric Levy models (GLMs). Here we consider the one- dimensional case. In what 
follows we shall assume that the reader is familiar with basic aspects of the theory of Levy 
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processes, as represented for example in Appelbaum (2004), Bertoin (2004), Cont & Tankov 
(2004), Kyprianou (2006), Protter (1990), Sato (1999), or Schoutens (2004). We recall that 
a Levy process on a probability space (Q, J 7 , P) is a process {X t } such that X = 0, X t — X s 
is independent of J^ s for t > s (independent increments), and 

F(X t -X s <y) = F(X t+h - X s+h < y) (4) 

(stationary increments). Here {Ft] denotes the augmented filtration generated by {X t }. 
For {X t } to give rise to a GLM, we require that 

E[e aXt ] < oo (5) 

for all t > 0, for a in some connected interval A containing the origin. Henceforth we 
consider Levy processes satisfying such a moment condition. It follows by the stationary 
and independent increments property that there exists a function ip(a), the so-called Levy 
exponent, such that 

E[e aXt ] = e t1p( - a) (6) 

for a G A. The process {M t } defined by 

M t = e aXt -^ (a) (7) 

is then called the geometric Levy martingale (or Esscher martingale) associated with {X t }, 
with volatility a. Indeed, by the stationary and independent increments property we have: 

E s [M t ] = E s [e aXt ] e~ ma) = E s [e a(Xt - Xa) } e QX *-*^ = e^)^ e <*x*-W<*) = m s . (8) 

Our geometric Levy model for asset prices will be put together as follows. First we 
construct the pricing kernel {vr f } f > . Let A > and assume that —A G A. Then set 

nt = e -rt e -XX t -t^X)^ (g) 

For a consistent pricing theory we require that the product of the pricing kernel and the 
asset price should be a martingale, which we shall assume is of the form 

n tSt = SoePX'-WW (10) 

for some /3 G A. From the formulae above we deduce that 

where a = (3 + A. We shall assume that a > and that a G A. It follows that the asset 
price can be expressed in the form 

S t = S e rt e R ^e° Xt -W a \ (12) 

where 

R(X, a) = z/>(cr) + ^(-A) - if>(cr - A). (13) 
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III. ON THE RISK PREMIUM ASSOCIATED WITH LEVY MODELS 

One sees that the function R(o~, A) is the risk premium, that is to say, the excess rate of 
return above the interest rate. Indeed, we have: 

E[3 t ] = 3 e rt+ ^y t . (14) 

The following result establishes a rather general property of geometric Levy models, and is 
indicative of why such models are of interest. 

Proposition 1. The excess rate of return in a geometric Levy model is positive, and is 
increasing with respect to the risk aversion and the volatility. 

Proof. We have ip( a ) = hi IE [e aXt ] , and thus 

1 E \{X t - X t ) 2 e aXt ] - E \X t e aXt ] 

ip (ol) — — lV * . t\ } -, where X t = J; . * . J . 15 

Formula f|T5|) shows that ip"{ot) > 0, and thus that the Levy exponent is convex as a function 
of a. Indeed, for any random variable £ satisfying m(a) := E [e a ^] < oo for a in some interval 
containing the origin, lnm(a) is convex (see, e.g., Billingsley 1995). Now consider four values 
of a in A such that a± < a 2 < a 3 < a 4, and for some h > suppose that 0:3 = a± + h and 
«4 = a 2 + h. Then we have: 

V>(«i) + -0(0:4) > ^(012) + ^(d 3 ). (16) 

To derive this inequality we note that the convexity of ip{ct) implies that ip ; (x+h) —tp'(x) > 
for x and x + h in A. Integrating ip'(x + h) — ip'(x) with respect to x between a\ and «2 one 
obtains fTTol) . Then since ^(0) = 0, and either —A <0<a — A<o"or —A < a — X < < a, 
it follows from fTTBT) by letting either h = a or h = A that 

0(-A)+0(o-) > V(o--A). (17) 

Therefore, -R(A, a) > 0. Furthermore, the convexity of the Levy exponent implies that 

9R ^ a) = jjj'(a - A) - t/j'(-X) > and dR ^ ^ = ^'(a) - t/j'(a - A) > 0. (18) 
o\ da 

Thus, R(\,o~) is increasing with respect to both A and a. □ 

We observe that the risk premium R(X, a) is in general a nonlinear function of the risk 
aversion (represented by A) and the risk (represented by a). This suggests that the notion of 
"market price of risk", so common in the finance literature, is somehow linked specifically to 
models based on Brownian motion, and is not quite the right idea in the context of general 
Levy models. Rather, risk premium is the more useful notion. In Section |V] we show that 
the only GLM leading to a bilinear risk premium is the GBM model. 

Properties of the risk premium can be examined further by use of the Levy-Khintchine 
representation for ijj(a>), which in the case of a Levy process admitting exponential moments 
takes the form 

1 r°° 

tp{a) = pa + -qa 2 + {e ax - 1 - ax\{\x\ < 1}) u(dx), (19) 
2 
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where p and q > are constants and u(dx) is the Levy measure (see, e.g., Sato 1999, 
Theorem 25.17). For any measureable set B G R, the expected rate at which jumps occur 
for which the jump size lies in the range B is v{B). It follows from ({T9~]) that the risk 
premium is given by 

/oo 
(e ax - 1)(1 -e- Ax >(dx), (20) 
-oo 

from which various of its properties can be deduced. In particular, the statement of Propo- 
sition [1] can be seen to follow from the fact that the argument of the integrand in ( 1201) 
is positive for a, A > 0, and that its first derivatives are positive. One can also calculate 
the higher derivatives of the risk premium with respect to risk aversion and volatility, and 
one deduces, for example, that in the case of a spectrally negative Levy process (downward 
jumps), we have d 2 a R < and d\R > 0, and that these inequalities are reversed in the 
case of a spectrally positive process. One also observes that, providing the tails of the Levy 
measure are not too fat, for small values of the risk aversion and volatility the risk premium 
is approximately bilinear. 

IV. LEVY MODELS FOR FOREIGN EXCHANGE 

When the geometric Levy model is extended to the case of foreign exchange, additional 
features arise that are of some interest. It is reasonable to require "numeraire symmetry" in 
the sense that if, for example, the dollar price of one pound sterling offers a positive excess 
rate of return above the interest rate differential, then the sterling price of one dollar should 
offer a positive excess rate of return above the reverse interest rate differential. 

We examine the GBM case first, where the situation is transparent. Let the dollar be 
the domestic currency, and the pound the foreign currency. Let St denote the price of one 
pound in dollars, and St the price of one dollar in pounds. We write r for the domestic 
(dollar) interest rate, and / for the foreign (sterling) interest rate, both assumed constant. 
Let A and a be positive constants, and let the dollar pricing kernel be given by ([T]). Then 
the GBM model for the foreign exchange rate (the dollar price of one pound) is given by 

S t = S e ir ~ f)t e^e^"^ 2 '. (21) 

We observe that the excess rate of return above the interest rate differential r — f is the 
product Aa, which is positive. For the corresponding inverse exchange rate we obtain 

S t = S e (/ - r)t e ( °- x)at e-^-i- 2 *. (22) 

In the case of St, we see that the excess rate of return above the reverse interest rate 
differential / — r is positive if and only if a > A. In equilibrium, we may presume that 
investors on both sides of the Atlantic wish to see the exchange rate promising a positive 
excess rate of return. The argument above shows that in a GBM model this possibility can 
be realised if a > A. This is the essence of the so-called Siegel (1972) paradox. 

Now let us look at the analogous situation in the context of a geometric Levy model. In 
the case of a GLM our elementary model for the exchange rate takes the form 



S t = S eC'-fl*e^ A ' ff >V**-^, 



(23) 
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where A and a are positive constants, and R(X, a) is given by (|T3|) . as before. Thus for the 
inverse exchange rate we obtain 

S t = S e (/ " r)t e Aix ' a)t Q-° x t-^(-°) , (24) 

where R(X, a) = —R(X, a) + ip(a) + ip(—a). A short calculation shows that 

R(\,a)=if>(-a)+il>(a-\)-iJ)(-\). (25) 

Now suppose that a > A. Then —a < — X < < a — A, which by taking ft, = a in equation 
(I16p implies that R(X, a) > 0. On the other hand, suppose that a < X. Then since either 
—A < —a < a — X < or —A < a — X < —a < 0, it follows from ffl6|) by letting either ft = cr 
or ft = A that -R(A, cr) < 0. One thus deduces the following: 

Proposition 2. // the volatility exceeds the risk aversion, then both (a) the excess rate of 
return on the FX rate, and (b) the excess rate of return on the inverse FX rate, are positive 
in a geometric Levy model for foreign exchange. 

One observes that the volatility of the inverse exchange rate in (I24p . as matters stand, 
appears to be negative, which is not consistent with our original characterisation of an 
investment-grade asset. We can however put the inverse exchange rate into "canonical" form 
by regarding it as being driven by the mirror process X t = —X t . Then a is the volatility, 
and we can regard the quantity A = a — X as being the associated risk aversion parameter 
for foreign agents. Writing b = cr, and regarding cr and A as independent variables, one can 
check that the inverse excess rate of return is monotonically increasing in both the volatility 
a and the foreign risk aversion factor A. Indeed, if we write ip(a) = i/)(—a) for the Levy 
exponent associated with the mirror process, then it is an exercise to verify that 

R(X,o-) = 4>(a) + i>(-X)-ij(a-X), (26) 

and one sees that the form of the foreign excess rate of return, when expressed in terms 
of the relevant Levy exponent, is identical to that of the domestic excess rate of return. 
The requirement a > X can be understood as an assertion that the foreign risk aversion is 
positive. Thus one might regard this constraint as a necessary feature of the model. 

Proposition 3. In a geometric Levy model for foreign exchange with a > X > , the excess 
rate of return on the inverse FX rate is increasing with respect to the independent variables 
A = cr — A and a = a . 

V. EXAMPLES OF GEOMETRIC LEVY MODELS 

It will be instructive to look at various explicit examples of geometric Levy models for asset 
prices, noting in particular the structure of the excess rate of return function in each case. 

Example 1: Brownian motion. In the case of a standard geometric Brownian motion 
model the Levy exponent is given by i^{a) = |a 2 , and hence 



R(X, a) = 4(a) + ip(-X) - ^(a - A) = l -a 2 + l -X 2 - l -{a - A) 2 = crA, (27) 
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which is positive. Further, for the inverse excess rate of return we have 

R(X, a) = M-a) + ifta - A) - i/j(-X) = 1 -a 2 + 1 -(a- A) 2 - l -\ 2 = a{p - A). (28) 

Thus R(X, a) is positive if and only if a > A. We observe that if one treats the quantities 
er = a and A = a — A as independent variables, then the inverse excess rate of return function 
is increasing with respect to each. 

One can ask to what extent the bilinear form of the excess rate of return determines 
the underlying Levy process. Is it uniquely Brownian motion that has this property? If we 
consider expression ffT3]) and set R(X, a) = Xa, then by taking two derivatives we deduce 
that tp"(a) = 1. Integrating twice and imposing the condition ^(0) = 0, we conclude that 
the general form of the Levy exponent compatible with a bilinear excess rate of return is 

ifj(a) =pa + ±a 2 , (29) 

where p is a constant, and we obtain a standard Brownian motion with drift. But the 
addition of a drift to the driving Levy process is irrelevant to the resulting pricing model in 
our scheme, since it cancels out in the formula for the geometric Levy martingale. Thus: 

Proposition 4. The geometric Brownian motion model is the only geometric Levy model 
with an excess rate of return that is bilinear in the risk aversion and the volatility. 

In fact, in any GLM the excess rate of return function is sufficient to determine the 
driving Levy process, up to an irrelevant drift. To establish that this is the case, we observe 
that if we differentiate each side of ( TT3|) with respect to A and a, and then take the limit as 
A approaches zero, we obtain 



(30) 

A=0 



dXda 

Integrating twice, and fixing the constant, we obtain the Levy exponent, modulo a drift. 

Example 2: Poisson process. Let {N t } be a standard Poisson process with jump rate 
m > 0. Then for any nonnegative integer n the distribution of N t is given by 

F{N t = n)= e'^^L. (31) 
n ! 

It follows that E[X(] = mt, and that the Levy exponent is ip(a) = m(e a — 1). The associated 
geometric Levy martingale with volatility a in this example is thus 

M t =exp{aN t -mt(e a -1)). (32) 

A calculation then shows that the excess rate of return function is manifestly positive, and 
increasing with respect to its arguments: 

R(X, a) = m(l — e _A )(e CT — 1). (33) 



We remark that since the jumps in the geometric Poisson model are upward, the "risk" that 
an investor faces is that there may be fewer jumps than one hopes for. This is made evident 
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if we combine the expressions for the geometric martingale and the excess rate of return 
function to obtain the following formula for the price of a non-dividend-paying asset: 

S t = S exp [rt + aN t - mte~ x (e a - 1)]. (34) 

Thus, the effect of investor risk aversion is to reduce the downward drift rate in the com- 
pensator term by attaching the factor e _A to it. For the associated pricing kernel one has 

ir t = exp[-rt- XN t -mt(e~ x -1)}, (35) 

and it is an exercise to check that the product of n t and St gives a geometric Poisson 
martingale with volatility o — A. In the event that St represents the price of a unit of foreign 
currency, then we include the foreign interest rate by setting 

S t = S exp [(r - f)t + aN t ~ mte~ x (e a - 1)]. (36) 

For the corresponding inverse exchange rate we obtain 

S t = S exp [(/ - r)t + R(X, a) - aN t + mt(l- e~ a )], (37) 

where R(X,a) = m(e cr ~ x — 1)(1 — e _<T ), in agreement with equation fl25l) . If a > A then the 
excess rate of return of the inverse exchange rate is evidently positive, and has the property 
of being increasing with respect to the independent variables a and a — A. 

Example 3: Compound Poisson process. Let {N t } be a standard Poisson process with 
rate m, and let {Yfc}fc € pj be a collection of identical independent copies of a random variable 
Y with the property that 

<f){a) := E [e aY ] < oo (38) 

for a in some connected interval A containing the origin. Writing 1{— } for the indicator 
function, one can check that 

oo 

X t = Y,Hk<N t }Y k (39) 
k=i 

defines a Levy process, and that the associated Levy exponent is given by the formula 
ij)(a) = m (<f)(a) — 1). In this example the excess rate of return function is 

R(X, a) = m {<p{a) + 0(-A) - <j>{a - A) - 1) , (40) 

and the fact that it is positive and is bi-monotonic in its arguments is evident as a conse- 
quence of the convexity of <f)(a), which follows from fl38|) . If St represents the price of a unit 
of foreign currency, then the resulting geometric Levy model for the exchange rate is 

S t = S exp [(r - f)t + aX t + mt (0(-A) - 0(a - A))]. (41) 

Example 4: Jump diffusion process. This example is a generalisation of the Merton 
(1976) jump diffusion model. The driver is a vector Levy process, one component being 
a standard Brownian motion, and the other a compound Poisson process with normally 
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distributed jumps. Write A and a for the risk aversion and volatility of the Brownian 
component, and (3 and 9 for the risk aversion and volatility of the jump component. Let Y 
be normally-distributed with mean zero and variance one. Then <p{a) = exp(|a 2 ), and the 
excess rate of return function is of the form 

R(X, a, f3,9) = Xa + m (e^ + e^ 2 - e^) 2 - l) , (42) 

which is positive and is monotonic in each of the variables. In Merton (1976) a key notion 
used to price options is the idea that there is no risk premium offered by the market in 
connection with idiosyncratic firm-specific risk. Merton assumes that jump risk is purely 
idiosyncratic and can be diversified away by holding a suitably broad portfolio. Merton 
argues that since the risk can be diversified away, the market awards no risk premium to 
investors who hold such assets. From a modern point of view, the assumption that jump 
risk is necessarily idiosyncratic is questionable: this is one of the lessons of the 2008 credit 
crisis. In our version of Merton's model, however, jump risk is being priced. 

Example 5: Gamma process. By a standard gamma process with growth rate m, we 
mean a process {7f} that has gamma-distributed stationary and independent increments, 
and satisfies E [7^ = mt and Var [7^ = mt. The density of 7 f is given by 

P( Tt E dx) = t{x > 0} X r[ J^ dx, (43) 

where T[a] denotes the gamma function. The identity V[a + 1] = aT[a] implies that the 
mean of % is thus justifying the interpretation of m as the growth rate. The associated 
moment generating function is 

E[e Q7t ] = (1 - a)- mt = e - m * ln d-«) ; (44) 
and hence the Levy exponent is 

ip(a) — -mln(l — a), (45) 

which is well-defined for a < 1. For some applications it is useful to consider the two- 
parameter family of so-called scaled gamma processes. By a scaled gamma process with 
growth rate \x and variance rate v 2 we mean a process {r t }o<t<oo with stationary and inde- 
pendent increments such that r = and such that T t has a gamma distribution with mean 
fit and variance u 2 t, where fi and v are parameters. Setting m = fJ 2 /^ 2 and k = v 2 / '//, one 
finds that n = Km and v 2 = n 2 m. One can thus think of m as a standardised growth rate, 
and scale parameter. The density of r t is given by 

¥(r t G di) = l{x > 0} ^ dx. (46) 

T[mt\ 

For fixed t the product mt is the so-called shape parameter of the distribution of the random 
variable r t . If {^ t } is a standard gamma process with growth rate m, then the process {r t } 
defined by r t = K^ t is evidently a scaled gamma process with standardised growth rate m 
and scale parameter k,, and for its moment generating function we have 

E[e aPt ] = (1 - aK)~ mt . (47) 
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See Dufresne et al. (1991), Dickson & Waters (1993), Madan et al. (1998), Baxter (2007), 
Yor (2007), and Brody et al. (2008) for various aspects of the theory of the gamma process. 

Now let {74} be a standard gamma process with growth rate m, and let a be a constant 
such that < a < 1. The associated geometric Levy martingale takes the form 

M t = (1 -a) ,T V 7t . (48) 

The jumps are upward, and the compensator is a deterministic decreasing process. If follows 
from ( JTBl and ( )45l) that the excess rate of return function is of the form 

If one takes the difference between the numerator and the denominator in the argument of 
the logarithm in (1491) . the result is (1 — a + A) — (1 — <r)(l + A) = aX , which is positive. It 
follows that R(X, a) > 0. By (|2"3"|) . the corresponding model for the foreign exchange rate is 

(\ mt 
1_ TTa) e<T7t ' (50) 

Variants of the geometric gamma model appear in Heston (1993), Gerber & Shiu (1994), 
and Chan (1999). We observe that the effect of risk aversion is to reduce the rate at which 
the compensator decreases, thus encouraging investors who might otherwise be concerned 
over the possibility of an insufficient rise in the underlying gamma process. For small A and 
u the risk premium is given approximately by Xa. 

As for the excess rate of return associated with the inverse exchange rate, by use of f[2"5]) 
we obtain the following expression: 

R(X,a) =mln- -. (51) 

We observe that (1 + A) — (l + er)(l — cr + A) = <j(o — A), from which it follows that R(X, a) > 
if and only if a > X. Numeraire symmetry thus imposes a bound on the risk aversion factor 
and we have < A < o < 1. Inverting (!50j) . and writing a = a and A = a — X, we find that 

,(f-r)t / i , ° \ --S-K 



S t = S e^~ r)t \1 + - — e- CT7t . (52) 

Thus St is driven by a negative gamma process, which jumps downward, and the compensator 
is a deterministic increasing process. Since < A < 1, the effect of foreign risk aversion is 
to increase the rate at which the compensator increases. For small values of A and a the 
inverse excess rate of return is given approximately by A a. 

Example 6: Variance gamma process. It will be convenient first to discuss the sym- 
metric variance gamma (VG) process. This is the process considered by Madan & Senata 
(1990) and Madan & Milne (1991). Then in the next example we discuss the more general 
asymmetric or "drifted" VG process of Madan et al. (1998). Both of these processes are of 
interest from a mathematical perspective and as a basis for financial modelling. There is a 
further extension of the model, due to Carr et al. (2002), which will not be discussed here. 
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The VG model relies on the use of a gamma process as a subordinator. Thus we begin with a 
standard gamma process {7 t } with rate m, and give it the dimensionality of time by dividing 
it by m. In this way we define a scaled gamma process {r t } by setting r t = m _1 7 t , and we 
observe that E [r t ] = t. We call {A} a standard gamma subordinator. The symmetric VG 
process {Vt}, with parameter m, is defined by letting {Wt} be a standard Brownian motion 
and setting Vt = Wr t - The associated moment generating function is thus 



E [exp(aVi)] = E [exp(aW r J] = E 



exp (^t 2 r t 



cr 
2m 



-mt 



(53) 



which is defined for a 2 < 2m. Clearly a must have units of inverse square-root time, since m 
has units of inverse time; but this is consistent with the fact that Vt has units of square-root 
time, like the Wiener process. The associated Levy exponent is 



ijt (a) 



a 



2 \ 



-m In 1 

2m J 



(54) 



and one can check the convexity of ip(a) in this example by observing that 



rj/'{a) = m(l 



a 

2m 



(55) 



As a consequence the geometric Levy martingale in the symmetric VG case takes the form 

(2 \ m t 
1 - |H exp (aW rt ) , (56) 

and the excess rate of return function, which is positive and monotonic, is 



R(X, a) — mln 



2m I \ 2m 



2m J 



(57) 



The corresponding VG foreign exchange rate is thus given by 



S t = S e^^ 1 



(a-Xf 
2m 



2 \ rni 



X 2 \ 



-mt 



exp(,W r ,). 



(58) 



We remark that in the case of the VG model one finds by use of fl55l) that the risk premium 
satisfies d 2 R > iff a > \a - A|, and d\R > iff A > \a - A|. 

A well known alternative characterisation of the VG process is as follows. Let {7^} 
and {7 t 2 } be a pair of independent standard gamma processes, each with rate m. Then 
the process defined by the difference between these two processes has both upward and 
downward jumps, and is symmetrical about the origin in distribution, with mean zero. If 
we normalise the difference by setting 



V 



2m 



W-7* 2 ) 



(59) 
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then it is easy to check that the variance of V t is t, and so we get a pure jump process that 
has some properties in common with Brownian motion. Indeed if we consider the moment 
generating function we find by virtue of the independence of the two gamma processes that 



E [exp(aV t )] 




-mt 



(60) 



and it is evident (Madan & Senata 1990) that (1591) has the law of a VG process. For large 
values of m, the distribution of Vt is approximately Gaussian. In particular, we have: 



( -^oi 2 t 



(61) 



Example 7: Asymmetric VG process. The representation of the VG process as the nor- 
malised difference between two independent gamma processes suggests two generalisations. 
One is that of Madan et al. (1998), where we consider an asymmetric difference between two 
independent standard gamma processes. Thus writing 

U t = k x1 ] - K 2l l (62) 

where K\ and k 2 are nonnegative constants, a calculation of the respective moment generating 
functions shows that U t is identical in law to a "drifted" VG process of the form 

U t = fir t + aW rt , (63) 

where /i and a are constants. The relationship between /x, a, k 1; k 2 , and m is given by 
fi = m(Ki — k 2 ) and a 2 = 2m/t 1 /t 2 , together with 

Ki = ( /i + a/u 2 + 2ma 2 \ and k 2 = ( — fi + a/u 2 + 2mff 2 ) . (64) 

2m V / 2m V / 

The Levy exponent ip{a) = — mln (1 — — n 2 ) at — K\K 2 a 2 ), which can be worked out by 
use of ( )47|) . can be equivalently written in the form 

ip(a) = —m In ( 1 — — a a 2 ] , (65) 

\ m 2m J 

where the range of a is — 1/k 2 < a < 1/ K\. It is straightforward to write down the associated 
excess rate of return function, and the corresponding expression for an asset price. In this 
example there is a single risk aversion factor. 

On the other hand, one can also envisage the situation where the two gamma drivers 
are regarded as separate sources of risk, each being assessed independently by the market. 
This situation arises in instances where investors are for some reason more worried about 
downward jumps than upward ones. More specifically, let us suppose that investors are 
more concerned about excessive losses than about insufficient gains. It is said that studies 
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in behavioural finance suggest that this may actually be the case. One can model such a 
situation rigorously by introducing an asymmetric pricing kernel of the form 

n t = e ~ rt (1 - AO" 1 ' (1 + A 2 ) m *e" Al7 'e A27 S (66) 

and an asset price process of the form 

/ \ mt / \ mt 

s - s ° e T-iTA7) ( 1 + r^) e °"'' e ^'- (67) 

Thus we have separate risk aversion factors for the upward jumps and the downward jumps. 
It is interesting to observe that in the case of "behavioural asymmetry" both the asset price 
and the pricing kernel are driven by extended VG processes — but there are two distinct such 
processes, one driving the pricing kernel, and the other driving the asset price. Indeed, the 
pricing kernel is driven by A17* — A27 2 , whereas the asset price is driven by 017* — 027^. 
These processes are synchronised in the sense that the times of their upward and downward 
jumps coincide and the magnitudes are proportional for a given jump type. Alternatively, 
we can model the two driving processes as different linear combinations of F t and Wp t - 

Example 8: Negative binomial process. This process is a species of compound Poisson 
process, and can be viewed as a special case of Example 3. It is nevertheless of considerable 
interest in its own right. See, for example, Kozubowski & Krzysztof (2009) for a general 
overview. The negative binomial process has a Levy exponent of the form 

^ a ) = mln (rr^)> (68) 

where < q < 1, and m > is a rate parameter. A short calculation shows that 



i?(A, a) = m In 



(1 - qe a )(l - qe~ 



(69) 



To see explicitly that this is positive we take the difference between the numerator and the 
denominator inside the logarithm to obtain q(e a — 1)(1 — e _A ), which clearly is positive since 
cr, A > 0. It follows that the argument of the logarithm is larger than one, and we have 
R(X, cr) > 0. For the inverse excess rate of return we obtain 



R(\, cr) = m In 



1 - q)(l - qe 



-X') 



(1 - qe~ a )(l - qe a ~ 



(70) 



To analyse the positivity of -R(A, cr) we again take the difference between the numerator and 
the denominator in the argument of the logarithm in flTOj) . We get c/(e°" _A — 1)(1 — e _<J ), 
which is positive only if a > A. 

The jumps of the negative binomial process are positive integers. There are two distinct 
representations for the process. The first of these takes the form of a compound Poisson 
process with the following characteristics: (a) the jump sizes have a so-called logarithmic 
distribution, given in the notation of Example 3 by 



(Y = n) = - / -q n ; 

ln(l — q) n 



(71) 
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and (b) the intensity fi of the underlying Poisson process is of the form \x = — mln(l — q). 
A straightforward calculation of the moment generating function of Y gives 

ln(l - qe a ) 

^ = ln(l - q) • (72) 

By the general theory of the compound Poisson process we know that ijj(a) = fi(<j)(a) — 1), 
which immediately leads to the Levy exponent fl68|) . If we write X t for the value of the 
negative binomial process at time t, we find that its probability mass function is given by 

which is the negative binomial distribution; and for the geometric Levy martingale we have 

M t=(^rrfy ^ aXt - (74) 

The second representation of the negative binomial process is reminiscent of the theory 
of the VG process. We take a standard Poisson process, with intensity 

A = mir A_, (75) 

where < q < 1 and m > as before, and subordinate it with a gamma subordinator {r t } 
with standardised rate parameter m. Thus the expectation of r t is t, and its variance is 
t/m, as in the theory of the VG process. The associated moment generating function is thus 

E [exp(aJVr t )] = E [exp(A(e« - 1)2})] = (l - ^=^) = (xf^)^ > ( 76 ) 

by virtue of the chosen intensity ( 1751) . and we are lead directly to the Levy exponent ( 1681) . 

It is natural to ask in the context of these various examples what information can be 
extracted (or "implied") about the values of model parameters when one is given option 
prices. In the case of the GBM model, for example, it is known that one can infer the 
value of volatility a, but that the option price is independent of risk aversion A. This can 
be checked directly by working out the price of a call option with strike K and expiry T 
by inserting ([I]) and (J2J) into the valuation formula Co = E [kt {St — K) + }. In a general 
geometric Levy model, this is no longer the case: option prices depend on both the risk 
aversion and the volatility. Indeed, a variety of different situations can arise, each with its 
own character. Thus in the Poisson model, there are two nontrivial model parameters — the 
risk aversion, and the jump rate m (the volatility is easily determined by observation of the 
price process); and a calculation shows that option prices depend on me _A , but not on m or 
A separately. Thus if we can estimate the value of the actual jump rate m by observations of 
the asset price, then A can be inferred from option prices. In the case of the gamma model, 
there are three nontrivial model parameters — the risk aversion, the volatility, and the jump 
rate. A calculation shows that option prices depend on m and on cr/(l + A), but not on a 
and A separately, so neither A nor a can be determined exactly from option prices. 
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VI. DIVIDEND PAYING ASSETS 



Thus far we have considered the case of non-dividend paying assets. From a conceptual 
point of view it is better, however, to think of an asset price as being determined by the 
dividend stream or cash flow produced by the asset. Hence with the inclusion of dividends 
the pricing model is characterised by: (a) a pricing kernel {vr t }; and (b) the dividend stream 
{D t } generated by the asset. The value of the asset at time t is regarded as a derived 
quantity which can be worked out by means of the fundamental relation 



S t =-E t 



oo 



(77) 



With this in mind, let us consider how one extends the GBM model when dividends are 
included. The answer to this problem is well known, but rather than assuming the conclusion 
we shall derive it from first principles by modelling the dividend stream and the pricing 
kernel, and working out the resulting price process for the asset. For a typical investment 
asset in the GBM situation we model the dividend stream by setting 

D t = D Q e^e am -^ H , (78) 

where Dq is the initial rate at which dividends are paid, 7 is the growth rate of the dividend, 
and the constant a > characterises the volatility of the dividend rate. We shall assume 
that the pricing kernel is of the form (JTJ. Substituting (JTJ and f!78|) into equation (ITT)) and 
performing a short calculation under the assumption that r + Act > 7, we deduce that 

S t = 1 D e 7*+<^-!<A (79 ) 

r + Act — 7 

Thus we obtain a stochastic generalisation of the Gordon (1959) growth model, and at time 
zero we have the following valuation formula: 

S = (80) 
r + Act — 7 

We observe that an increase in risk aversion has the effect of lowering the asset price, 
everything else being the same. Defining the proportional dividend rate by setting 5 = 
r + Act — 7, which by assumption is positive, we are able to deduce the relation D t = SS t , 
and we find that the asset price process is given by 

S t = S e ir - s+X(T)t e aWt -^ 2 \ (81) 

As we stated earlier, the resulting expressions for {St} and {D t } are of course familiar: the 
point is that we derive these formulae here rather than assume them. 

In the situation of a geometric Levy model it is remarkable that essentially the same line 
of argument carries through. Thus we assume a pricing kernel of the form (j5J), and a volatile 
dividend stream of the form 

D t = D e 7 V*'-^ (<7) . (82) 
It is then an exercise to check that the fundamental relation gives 

S t = So e {r - 5)t e R{x ' a)t e'X'-W") , (83) 



16 



and a proportional dividend flow D t = 5S t , where 5 = r + R(\, o) — 7- The resulting initial 
valuation formula is 

s ° = ; + R °U -7 • (84) 

This relation ties together in the context of a general geometric Levy model the values of 
the initial asset price, the initial dividend rate, the interest rate, the risk aversion level, the 
dividend volatility, and the dividend growth rate. 



VII. MULTI-FACTOR MODELS WITH PREDICTABLE VOLATILITY 

There are several reasons for extending the analysis to higher dimensions. First, we would 
like to consider models for a market consisting of a number of different assets. Second, even 
in the consideration of a single asset, it is natural to introduce the additional complexity 
of a higher-dimensional process to describe its dynamics. Both situations are familiar in 
the context of Brownian motion driven models. It seems to be advantageous to envisage 
the entire market as being driven by a single higher-dimensional Levy process. We can use 
essentially the same notation as in the one-dimensional case. Now {X t } is understood to be 
a vector Levy process. For the Levy exponent we still have (161) , but now a is understood 
to be a vector, and there is an implicit inner product between a and X t in the exponent on 
the left hand side of ([6]). The Levy exponent is a function of the n components of a. 

The model is thus determined as in Section HJ with the assumption that the market 
filtration is generated by a vector Levy process, and with the specification of (i) a pricing 
kernel, and (ii) a collection of investment grade assets, driven collectively by {X t }. We are 
accustomed, in the mult i- dimensional Brownian case, to regard such a higher dimensional 
driver as being built from a set of independent drivers that can be isolated after a suitable 
linear transformation. No such simplification is readily at hand for a general vector Levy 
process. Nevertheless for applications it is useful to consider the case where the components 
are assumed to be independent. This encompasses a large class of models, including the 
higher- dimensional Brownian motion models. 

The pricing kernel is a process of the form (Q where A is now understood to be a vector 
risk aversion factor. If the components of the vector Levy process are independent, the 
Levy exponent separates into a sum of terms, one for each component of its argument, each 
term being the marginal Levy exponent associated with one of the risk factors. Next we 
introduce a set of investment-grade assets, each of the form 012p for some choice of the 
vector volatility a. We require that A and a are "positive" vectors — that is, they belong to 
the cone of vectors with the property that all components are nonnegative and at least one 
component is positive. For a generic asset the excess rate of return takes the form (1131) . only 
now the arguments are understood to be vectors. When the Levy process has independent 
components, the excess rate of return separates into a sum of terms, each being the excess 
rate of return associated with one of the components. In that case, we see by use of the 
arguments presented earlier that each term is nonnegative, and at least one is positive; as 
a consequence the total excess rate of return function is positive. Similarly one sees that 
in the case of independent components the excess rate of return is increasing with respect 
to the individual components of the risk aversion vector and the volatility vector. Thus we 
obtain the following result. 
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Proposition 5. The excess rate of return function R(X, o~) in a multi- dimensional geometric 
Levy model with independent Levy drivers is positive, and is increasing with respect to each 
of the components of the risk aversion vector X and the volatility vector a. 

In the discussion so far, we have assumed for simplicity that the interest rate, the risk 
aversion, and the asset price volatilities are constant. Indeed, as with many financial models, 
various characteristic features of the model are already present under the assumption of 
constant coefficients; but for practical applications, and to take the theory further, we need 
to relax this condition. Thus in the case of geometric Brownian motion models we consider 
the situation where the risk aversion {X t } and volatility {a t } are adapted vector-valued 
processes, and are chosen in such a way as to ensure that the process defined by the expression 

M t = exp (J a s dW s - \ J a 2 s ds^j (85) 

is an {J^j-martingale for a t = —X t , a t = o~ t , and a t = a t — X t ■ It suffices that {a t } should 
be bounded. More generally, we consider the situation where {Xt} is a Levy process with 
exponential moments, {Ft} is the associated filtration, and {ct t } is a predictable process, 
adapted to {Ft}, chosen in such a way that a t G A for t > and that the local martingale 
defined by 

M t = exp a s dX s - jf ^(a.)ds^ (86) 

is a martingale. If a predictable vector process {at} satisfies these conditions then we say it 
is admissible. Thus we consider a market model of the following form. Let the exogenously 
specified short rate process {r t } be adapted to {Ft}, and be such that the unit-initialised 
money market account 



B t = exp\^J r s dsj (87) 

is finite almost surely for t > 0. Let the {J-jj-adapted vector risk aversion and volatility 
processes {X t } and {a t } be positive, and be such that the processes {—X t }, {o~t}, and {o- t — X t } 
are admissible in the sense described above. The pricing kernel is taken to be of the form 

7r t = exp (- 1^ r s ds -J X s dX s - ip(-X s ) ds^j . (88) 

The corresponding expression for the price of a typical non-dividend-paying asset is then 

S t — So exp I / r s ds+ R(X s ,a s )ds+ / a s dX s - / ip(a s ) ds ) , (89) 
\Jo Jo Jo Jo / 

where R(X, a) is the excess rate of return function associated with the given Levy exponent. 
Clearly both {irtB t } and {fttSt} are martingales. It should also be evident that the following 
statement holds, which is to be understood as an expression of the fact that the asset offers 
a rate of return greater than the interest rate. 
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Proposition 6. Let the asset price in a model with predictable volatility and risk aversion, 
driven by a vector of independent Levy processes, be given by (189|) . and define the money 
market account by flH7|) . Then the asset price, expressed in units of the money market 
account, is a submartingale. 

In the literature, stochastic volatility models are often introduced by the method of a time 
change. In the context of simple parametric volatility this amounts to the observation that if 
{W t } is a standard Brownian motion, then the processes defined by Z t = aW t and Z' t = W a i t 
have the same law. Thus instead of introducing the volatility as a coefficient measuring the 
sensitivity of the asset price to the underlying Brownian motion, it is introduced by "speeding 
up"( or "slowing down") the Brownian motion: the effect is equivalent. In the case of jump 
processes, the two transformations are clearly inequivalent; thus, for example, in the case 
of the Poisson process the effect of scaling the process (magnifying the the jumps) is quite 
different from that of scaling the time (speeding up the arrival of jumps). It seems that in 
the general situation one wishes to consider both effects. Thus if the market is driven by a 
vector Levy process {X t }, then we introduce a vector of sensitivity parameters a, as well as 
a time dilation factor c, and let a typical asset be driven by the process {aX ct }. Again, we 
think of the entire market as being driven by a single vector Levy process so it is consistent 
that there is but a single overall time change for that process. On the other hand if the 
time change is random then we introduce a subordinator {c t } and the time-changed process 
is given by {aX Ct }. If the subordinator is itself a Levy process then one stays within the 
category of models already under consideration — thus the VG process can be obtained by 
subordinating Brownian motion with a gamma process, and the negative binomial process 
can be obtained by subordinating a Poisson process with a gamma process, as we have seen. 

VIII. CONCLUDING REMARKS 

Geometric Levy models have a surprisingly wide range of desirable properties. As we have 
seen, once suitable inequalities are imposed on the volatility and risk aversion parameters, 
the convexity of the Levy exponent ensures that the excess rate of return function is positive 
and is monotonic. In foreign exchange models, numeraire symmetry can be ensured by 
imposing a further inequality on the relation between the volatility and the risk aversion. 
In the extended version of the model, the market is driven by a vector of independent 
Levy processes, and the risk aversion and volatility coefficients are taken to be predictable 
processes satisfying suitable integrability conditions, these conclusions remain valid. Our 
approach is based on use of the physical measure P. We emphasise the importance of the 
pricing kernel method, because this leads to a unified view of the role of Levy models in 
finance, allowing one to separate pricing issues from hedging issues. In particular, we make 
no use of the idea of "trying to find an equivalent martingale measure" by some recipe when 
one is given a set of price processes. Rather, the pricing kernel is to be regarded as an 
essential component of the theory from the beginning. One needs the pricing kernel before 
one can speak of prices, because the value of a security is determined by the random cash 
flows that it produces, and these need to be valued by use of the pricing kernel. In this respect 
our point of view diverges in spirit from the earlier literature on Levy models in finance, as 
represented by Gerber & Shiu (1994), Eberlein & Keller (1995), Eberlein & Jacod (1997), 
Chan (1999), Raible (2000), Kallsen & Shiryaev (2002), Fujiwara & Miyahara (2003), Esche 
&: Schweizer (2005), Eberlein et al. (2005), Hubalek & Sgarra (2006), and others. Exceptions 
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include: Madan & Milne (1991) who are able to identify the rate of return in their study of 
the VG model; Heston (1993) who with the introduction of a "pricing operator" in his study 
of the gamma model offers a point of view similar in some respects to ours; and Madan 
(2006) where risks are priced by such kernels in a particular equilibrium. 

The approach we have outlined for asset pricing in a Levy setting with predictable risk 
aversion and volatility is also useful in the theory of interest rates, since it allows one to 
generalise the HJM (Heath et al. 1992) framework in a natural way to the Levy category, 
without the need of introducing instantaneous forward rates, but in a way that guarantees 
positive excess rates of returns on bonds, and is formulated in the P-measure, making it 
suitable as a practical basis for risk management, forecasting, and scenario analysis. 
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